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TITLE OF THE INVENTION 
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This application claims priority of U.S. Provisional 
v Patent Application No. 60/476,847 filed June 9, 2003 
'entitled ESTIMATION OF BACKGROUND NOISE AND ITS EFFECT ON 
SONAR RANGE ESTIMATION, and U.S. Provisional p Patent 
15 Application No. 60/476,849 filed June 9, 2003 entitled 'NOISE 
ADAPTIVE SONAR SIGNAL PROCESSOR. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 

20 The present invention was made with Government support 

under U.S. Government Contract Nos. ARO DAAD 19-02-1-0403 
and ONR N00012-02-C-02960. 

BACKGROUND OF THE INVENTION 

25 The present application relates generally to sonar 

systems, and more specifically to systems and methods of 
increasing the accuracy of sonar range measurements. 

Sonar systems are known that employ sonar pulses 
reflected from an object or target to estimate a distance to 

30 the target (also known as estimating the range of the 
target) . Such sonar systems include side-scan sonar 
systems, forward and downward- loo king sonar systems, 2-D and 
3-D sonar systems, and synthetic aperture sonar systems. A 
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conventional sonar system for performing sonar range 
estimation is typically configured to transmit one or more 
sonar pulses comprising sonic or supersonic pressure waves 
toward a selected target , and to receive one or more sonar 
5 pulses reflected from the target. Such reflected sonar 
pulses, which are commonly called echoes or returns, may 
include a m Significant amount of background noise, and/or 
other interfering signals in addition to a reflected sonar 
Signal of interest. The conventional sonar system typically 
10 includes a coherent receiver (also known as a matched filter 
receiver) configured to receive both the echo and a 
representation of the transmitted sonar pulse. For. example, 
the coherent receiver may comprise a cross correlator. The 
echo and the representation of the transmitted sonar pulse 
15 are cross-correlated within the coherent receiver to 
generate a peak cross correlation value, which is compared 
to a predetermined threshold value. If the cross 

correlation value is greater than the predetermined 
threshold value, then the reflected sonar signal of interest 
20 has been successfully detected. The conventional sonar 
system then utilizes the cross correlation peak to obtain a 
. measure of the range of the target. 

One drawback of the above-described conventional sonar 
system is that the level of background noise and/or other 
25 interfering signals contained within the echo or return may 
be sufficient to cause the reflected sonar signal to go 
undetected or to be falsely detected, thereby causing the 
cross correlator to produce inaccurate range measurements. 
Such inaccurate range measurements are likely to occur in 
30 low signal-to-noise ratio (SNR) sonar environments, in which 
the noise power within the echo may be comparable to or 
greater than the reflected signal power. This can be 
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problematic in sonar range estimation systems because a 
reduction in the measurement accuracy of the cross 
correlator typically leads to a concomitant reduction in 
sonar range accuracy. 
5 Prior attempts to increase the accuracy of sonar range 

measurements in noisy sonar environments have included 
filtering out at least some of the background noise before 
providing the echo to the cross correlator. However,, such 
attempts have generally 'not worked well enough to allow 

10 successful detection of reflected sonar signals and accurate 
estimation of range in low SNR sonar environments. This is 
due, at least in* part, to the fact that sonar systems 
typically receive sonar pulses that include various types of 
noise from a variety of different noise sources. For 

15 example, a sonar system may transmit sonar pulses through a 
medium such as water from a ship or submarine that produces 
noise across a wide range of frequency. Further, other 
ships, submarines, or structures producing noise across wide 
frequency ranges may be within the vicinity of the sonar 

20 system. Moreover, the natural interaction of the water and 
objects within the water including the selected target may 
produce a substantial amount of ambient noise. 

Various adaptive beam forming and noise cancellation 
techniques have been employed in an attempt to minimize the 

25 effect of noise on sonar range estimation. However, such 
techniques generally require some knowledge of the 
characteristics of the noise source. For this reason, such 
techniques are usually, most effective when used to minimize 
the effect of the noise produced by the ship or submarine 

30 carrying the sonar system, which is relatively easy to 
characterize. However, such techniques are typically 
ineffective when the source of the noise is unknown. 
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In addition, sonar ranging systems may receive echoes 
from a plurality of selected (and unselected) targets, each 
target having its own associated noise level, and it may be 
desirable to determine the noise level and range of each 
5 target separately. Such noise associated with multiple 
targets may be stationary or non-stationary, linear or 
nonlinear, or additive or non-additive. Further, \ at least 
some of the background noise may result from reverberations 
andy'or random signal distortions of , the transmitted* .or 
10 reflected sonar pulse, and therefore both the noise level 
and its structure may be significantly affected by the 
'transmitted sonar signal. However, conventional sonar 
systems are generally incapable of accurately estimating, 
noise levels and target ranges in the presence of non- 
15 stationary, nonlinear, non-additive, and/or signal-dependent 
noise. 

Moreover, the density and temperature of the 
transmission medium (e.g., water) and the frequency of the 
transmitted/reflected sonar signal may affect the decay rate 

20 of the sonar pulse propagating through the medium. In 
addition, the absorption of certain frequencies of the 
transmitted sonar pulse by the target may affect the 
strength of the resulting echo or return. The sonar pulse 
strength, frequency range, medium dispersion properties, 

25 shape, and duration may also affect the accuracy of sonar 
range measurements. However, conventional sonar systems are 
generally incapable of fully compensating for such factors 
when called upon to generate accurate noise and range 
estimates . 

30 . It would therefore be desirable to have a system and 

method of performing sonar range estimation that takes into 
account the effects of background noise whether the noise is 
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stationary or non-stationary, linear or non-linear, additive 
or non-additive, or signal-dependent or non-signal- 
dependent. It would also be desirable to have a method of 
performing sonar range estimation that can be used to 
5 increase sonar ranging accuracy in low SNR sonar 
environments. 

i 

, BRIEF SUMMARY. OF THE INVENTION 

In accordance with the present invention, a system and 

10 method is provided that permits sonar range estimation in a 
noisy sonar environment. Benefits of the presently 

disclosed system and method are achieved by obtaining an 
estimate of the signal-to —noise ratio (SNR) in the sonar 
environment, determining the effect of the estimated SNR on 

15 sonar ranging accuracy, and adapting the system and method 
to account for the determined effect of the SNR. 

In one embodiment, a system for performing sonar range 
estimation in a noisy sonar environment comprises a sensor, 
a transmitter, a receiver, a plurality of band-pass filters, 

20 a cross correlator, and a data analyzer. The sonar 
transmitter is configured to transmit at least one sonar 
pulse through a transmission medium such as air, water, 
soil, living tissue, or any other suitable transmission 
medium. The transmitted pulse travels through the 

25 transmission medium until it strikes an object, which 
returns at least one reflected pulse (echo or return) to the 
sensor. The sensor is configured to provide a signal 
representative of the echo to the sonar receiver, which 
subsequently provides an indication of the echo to the 

30 plurality of band-pass filters. Each band-pass filter is 
configured to pass a respective sub-band of frequencies , in 
which the respective frequency sub-bands substantially span 
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the frequency range of the transmitted pulse. The echo and 
a representation of the transmitted pulse are filtered by 
the respective band-pass filters, which provide filtered 
versions of the echo and pulse to the cross correlator. The 
5 cross correlator . is configured to. perform multiple cross 
correlation operations, in which each cross correlation 
operation operates on filtered versions of the echo and 
pulse produced by a respective one of the band-pass filters. 
The cross correlator provides cross correlation output data 

10 to the data analyzer, which is operative to analyze the data 
to determine the variability of cross correlation peaks 
within eaqh frequency sub-band, and to" identify the lowest 
frequency sub-band having a corresponding low peak 
ambiguity. By referencing the peak variability versus SNR 

15 corresponding to the identified frequency sub-band, which 
may be theoretically or empirically determined, the data 
analyzer is further operative to provide an estimation of 
the SNR within the sonar environment of interest, and an 
indication of the center frequency of sonar pulses that may 

20 be employed to perform accurate sonar ranging within that 
environment. The sonar system adapts to the estimated SNR 
by dynamically changing the frequency of the transmitted 
pulses to concur with the center frequency indication 
provided by the data analyzer. 

25 By adapting to the estimated SNR within the sonar 

environment of interest, the presently disclosed system and 
method may be employed to obtain sonar range estimations 
within that environment with increased accuracy. 

The disclosed system and method provides an estimation 

30 of the instantaneous SNR of a signal (e.g., a sonar pulse or 
"ping" ). Specifically, if the ping strikes several layers 
inside the object, the return signal from the inner layers 
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is weaker and therefore the SNR is reduced. Due to 
reverberations inside the object, the SNR may be. further 
reduced. The presently disclosed system and method estimate 
the SNR for each returning ping separately, and therefore 
5 adapt to the correct bandwidth for each returning ping. 

Other features, functions, and aspects of the invention 
will be evident from the Detailed Description of the 
Invention that follows, 

■ « • 

.* «s 

10 " BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

The invention will be moire fully understood with 
reference to the following Detailed Description of the 
Invention in conjunction with the drawings of which: 

Fig. 1 is a block diagram of a system for performing 
15 sonar range estimation in a noisy sonar environment 
according to the present invention; 

Figs. 2a-2b are diagrams of ambiguity functions 
illustrating the effect of noise level on the variability of 
cross correlation peaks; 
20 Fig. 3a is a diagram illustrating peak variability as a 

function of signal-to-noise ratio and center frequency for a 
plurality of frequency sub-bands; 

Fig. 3b is a diagram illustrating a performance curve 
derived from the diagram of Fig. 3b; 
25 Figs. 4a-4c are diagrams illustrating empirical 

estimations of sonar ranging accuracy; and 

Figs. 5a-5b are flow diagrams illustrating methods of 
operating the system of Fig. 1. 



30 



DETAILED DESCRIPTION OF THE INVENTION 
U.S. Provisional Patent Application No. 60/47 6,847 
filed June 9, 2003 entitled ESTIMATION OF BACKGROUND NOISE 
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AND ITS EFFECT ON SONAR RANGE ESTIMATION, and U.S. 
Provisional Patent Application No. 60/476,849 filed June 9, 
2003 entitled NOISE ADAPTIVE SONAR SIGNAL PROCESSOR are 
incorporated herein by reference. 
5 Systems and methods of performing sonar range 

estimation are disclosed that may be employed to increase 
the accuracy of sonar ranging in low SNR sonar environments. 
The presently disclosed systems and methods are operative to 
estimate the signal-to-noise (SNR) within the sonar 
10 environment of interest, and to adapt to the estimated SNR 
to obtain sonar. range estimations with increased accuracy. 

Fig. 1 depicts an illustrative embodiment of a system 
100 in accordance with the present invention that permits 
sonar range estimation in a noisy sonar environment. In the 
15 illustrated embodiment, the system 100 comprises a sensor 
102, a transmitter 103, a receiver 104 , a plurality of band- 
pass filters 106, a cross correlator 108, a data analyzer 
110, and a controller 111. It is noted that the 

illustrative embodiment of the system 100 described herein 
20 is suitable for performing range estimation in a noisy sonar 
environment. For example, the sonar system 100 may be 
adapted for (1) marine exploration in an underwater 
environment, (2) seismic exploration in a soil environment, 
(3) medical ultrasound in an environment comprising living 
25 tissue, or any other suitable use in a sonar environment. 
It should be understood, however, that the presently 
disclosed system 100 for performing range estimation may 
also be adapted for use in radar systems, microwave systems, 
laser systems, or any other suitable system. 
30 Specifically, the sonar system 100 includes the sonar 

transmitter 103, which is configured to transmit at least 
one sonar pulse through a transmission medium such as water. 
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The transmitted sonar pulse travels through the water until 
it strikes an object or target 112 in the water, which 
returns at least one reflected sonar pulse (commonly known 
as an echo or return) toward the sonar sensor 102 . For 
5 example, the sonar sensor 102 may comprise one or more 
hydrophone sensors. The sensor 102 is configured to sense 
. the echo, and to provide a signal representative of the echo 
.to the sonar receiver 104, which in turn provides an 
indication of the echo to the plurality of band-pass filters 
10 106. 

In the illustrated embodiment, each of the band-pass 
filters 106. 1-106. n is configured to* pass a respective sub- 
band of frequencies, in which the respective frequency sub- 
bands are defined to span the frequency range of the sonar 
15 pulse transmitted by the sonar transmitter 103. For 
example, in the event the frequency range of the transmitted 
sonar pulse is about 87 kHz, the respective frequency sub- 
bands may approximately range from 6-18 kHz, 18-31 kHz, 31- 
43 kHz, 43-56 kHz, 56-68 kHz, 68-81 kHz, and 81-93 kHz, 
20 thereby spanning the pulse frequency range of 87 kHz. 
Accordingly, in this illustrative example, the center 
frequency of the frequency sub-bands may be approximately 
equal to 12 kHz, 25 kHz, 37 kHz, 50 kHz, 62 kHz, 75 kHz, and 
87 kHz, respectively. 
25 As shown in Fig. 1, the plurality of band-pass filters 

106 is configured to receive the echo indication from the 
sonar receiver 104, and to receive a representation of the 
sonar pulse transmitted by the sonar transmitter 103 over a 
line 120. The echo and the representation of the 
30 transmitted pulse are filtered by each of the respective 
band-pass filters 106. 1-106. n, which subsequently provide 
band-pass filtered versions of the echo and the transmitted 
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pulse to the cross correlator 108. The cross correlator 108 
is configured to perform multiple cross correlation 
operations on the filtered echoes and pulses. Specif ically, 
the cross correlator 108 cross-correlates the filtered 
5 . versions of the echoes and pulses provided by each of the 
band-pass filters 106. 1-106. n, and provides corresponding 
cross, correlation output data to the data analyzer. 110, 
which is operative to analyze the data to ' determine, the 
variability of cross correlation peaks within each frequency 

10 sub-band , and to identify the lowest frequency sub-band 
having a corresponding low peak ambiguity, thereby allowing 
the data analyzer 110 to estimate the SNR of the sonar 
environment. The operation of the data analyzer within the 
presently disclosed system is described in greater detail 

15 below. The data analyzer 110 is further, operative to 
provide information to the controller 111, which is 
operative to cause the sonar transmitter 103 to modify the. 
characteristics of the transmitted pulses to adapt to the 
estimated SNR. In this way, more accurate sonar ranging is 

20 achieved. 

The cross correlation of an echo and a pulse may be 
expressed as 

M>e o \J/p(T) - J\l/e(t)l|/p(t+T)dt => 

25 = J\|/p(t)\j/p(t+T+T 0 )dt + frfp(t)r|(t+T)dt, (1) 

in which the first term * J\j/ P (t) \|/ p (t+T+x 0 ) dt" is the auto- 
correlation of the pulse centered at time x 0 , the second 
term xx Jv|/ P (t)T) (t+x) dt" is representative of band-limited white 
30 noise with frequency limits defined by the spectrum of the 
pulse, and the integration operation in each term is 
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performed from -oo to +00. Figs. 2a-2b depict representative 
ambiguity functions 201-203 that may be employed to describe 
the output provided by the cross correlator 108 (see Fig. 
1) . Because the cross correlator 108 cross correlates the 
5 filtered echo and pulse provided by each of the band-pass 
filters 106. 1-106*. n, each of which passes a respective sub-. 
* band of frequencies, it* is understood that an' ambiguity, 
function may be constructed corresponding to each ' of the 
frequency sub-bands. 
10 As shown in Figs. 2a-2b, ^ the ambiguity functions 201- 

203 are expressed as functions of pulse amplitude (vertical 
axis, dB) and delay time (horizontal axis, seconds), which 
is proportional to range. Specifically,, the ambiguity 
functions 201-203 correspond to the cross correlation of 
15 respective echo and pulse pairs having approximately the 
same frequency range but different center frequency fc 
(i.e., mean integrated frequency). For example, the 
ambiguity function 201 corresponds to the cross correlation 
of a first echo and pulse pair having a low center frequency 
20 fcl, the ambiguity function 202 corresponds to the cross 
correlation of a second echo and pulse pair having an 
intermediate center frequency fc2, and the ambiguity 
function 203 corresponds to the cross correlation of a third 
echo and pulse pair having a high center frequency fc3. 
25 Fig. 2a depicts a detailed view of the main lobes of the 
ambiguity functions 201-203, and Fig. 2b depicts the main 
lobes and the side lobes of the ambiguity functions 201-203. 
Each one of the ambiguity functions 201-203 comprises a 
respective peak value, which is indicative of the range of 
30 the object or target returning the echo. 

In high SNR sonar environments (i.e., when the noise 
level is low) , the peak of the ambiguity function is 



WO 2005/006009 



PCT/US2004/018059 



generally located at the main lobe of the function. In this 
case, the peaks of the ambiguity functions 201-203 are 
regarded as having low ambiguity, and may be located within 
the width of the main lobes of the functions at about time 
5 x 0/ as. illustrated by the vertical line of Fig. 2a. It is 
appreciated that* the time x 0 corresponds to the actual range 
"of the target. The effect of the low level of , noise . in such 
high SNR sonar environments is to jitter the position of the 
peak around' the time t 0 . To a first approximation, the 
10 magnitude of this jitter (also known as peak variability) is 
relatively low, e.g., the peak variability is typically less 
than the width of the* main lobes,, as. illustrated by the 
horizontal lines of Fig. 2a. The lengths of the horizontal 
lines of Fig. 2a are indicative of the levels of peak 
15 variability associated with the respective ambiguity 
functions 201-203. In the illustrated embodiment, the 
lowest peak variability is associated with the ambiguity 
function 203 (high center frequency fc3) , and the highest 
peak variability is associated with the ambiguity function 
20 201 (low center frequency fcl) . 

In low SNR sonar environments (i.e., when the noise 
level is high, for example, when the noise level is of the 
order of the difference between the amplitudes of the main 
lobe and the first side lobe) , the peak of the ambiguity 
25 function may not be located within the main lobe of the 
function, but instead may be located at one of the side 
lobes. In this case, the peaks of the ambiguity functions 
201-203 are regarded as having high ambiguity, and may be 
located (1) within the width of a side lobe at about time x. 2 
30 for function 203, (2) within the width of a side lobe at 
about time x- 2 . 5 for function 202, and (3) within the width 
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of a side lobe at about time x_ 3 for function 201, as 
illustrated by the vertical lines 203a, 202a, and 201a, 
respectively, of Fig. 2b. The effect of the high level of 
noise in such low SNR sonar environments is to significantly 
5 increase the . peak variability, thereby increasing the 
.potential error in" sonar range estimation. The horizontal 
lines in Fig. 2b illustrate the potential error in range 
estimation^that can result from such high noise levels. 
r Fig.' 3a depicts peak variability as a function of SNR 
10 (dB) and* center frequency fc for a plurality of frequency 
• sub-bands. In this illustrative example, peak variability 
is expressed in terms of * root mean square error (RMSE, - 
seconds), which is a temporal representation of the 
potential error in range estimation. Further, the center 
15 frequencies fc of the frequency sub-bands are equal to 12 
kHz, 25 kHz, 37 kHz, 50 kHz, 62 kHz, 75 kHz, and 87 kHz, 
respectively, and the centralized root mean square bandwidth 
Bcrms of transmitted pulses is fixed at 2 . 1 kHz. 

For example, peak variability curves 301-307, as 
20 depicted in Fig. 3a, may be obtained via Monte Carlo 
simulations. Specifically, the transmitted pulses may be 
expressed as cosine packets of the form 



\|/«,t,(t) = K a , T ,exp(-t 2 /2a 2 )cos(27rnt) , (2) 

25 

in which w Tf is the center frequency; w a" is the standard 
deviation of a peak location in time, which controls the 
spread in time of the pulse and its frequency bandwidth; 
and, n Ko, n " is a normalization factor such that 

30 

Va,ti(t)dt = 1, (3) 
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in which the integration operation is performed from -co to 
+00. Further , white noise may be added to the pulse to 
generate a noisy echo for the simulation, and a temporal 
5 indication of the range estimate may be computed as the time 
.corresponding to the maximum amplitude of the cross 
correlation between the echo and pulse. 

As shown in Fig. 3a, each of the simulation curves 301- 
307 is approximately linear within a first SNR range of 
10 about 35-50 dB (see also region I of Fig. 3b) . Further , for 
each curve 301-307, there is a sharp transition from lower 
RMSE- levels . to higher RMSE levels within a second SNR range 
of about 15-35 dB (see also region II of Fig. 3b), thereby 
indicating significant increases in peak variability. 
15 Within a third SNR range of about 5-15 dB (see also region 
HI of Fig. 3b), the curves 301-307 are again approximately 
linear. It is noted that the curve 308 depicted in Fig. 3b 
is a performance curve comprising a partial composite of the 
peak variability curves 301-307, including break points 1-9 
20 (see Fig. 3a). Accordingly, as the SNR decreases (i.e., as 
the noise level increases), the corresponding RMSE values 
gradually increase within region I until sharp transitions 
occur from lower RMSE levels to significantly higher RMSE 
levels within region II - the RMSE values then continue to 
25 increase more rapidly at the higher RMSE levels within 
region III. It is noted that within region III, the sonar 
range resolution falls sharply to a point where the sonar is 
ineffective and the target is considered to be out-of -range . 
Specifically, within region I, the simulation curves 
30 301-307 approximately track a line 310 (see Fig. 3b), which 
may be defined as 
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a - ^TcBRMsdr 1 , (4) 

in which >x a" is the standard deviation of a peak location in 
time and is proportional to the RMSE, "Brms" is the root mean 
square bandwidth of the pulse, and "d" is the SNR. A 
derivation of . equation (4) is described in Probability and 
Information Theory . with Applications to Radar, P.M. 
Woodward, New York, McGraw-Hill Book Company, Inc., 
copyright 1953, which is "incorporated herein by reference. 
It is noted that Brms may be expressed as 

Bbms - (ff 2 Ps D (f)df) 1/2 , (5) 

in which "P S D(f)" is the power spectral density of the 
15 pulse, and the integration operation is performed from 0 to 
+oo. Further, d may be expressed as 

d = (2E/N 0 ) 1/2 , (6) 

20 in which "E" is the total energy of the echo, and XX N 0 " is 
the spectral density of the noise. Accordingly, 

SNR(dB) = 201ogi 0 d. (7) 

25 Moreover, following the sharp transitions from lower 

RMSE levels to higher RMSE levels within region II (see Fig. 
3b) , the simulation curves 301-307 approximately track a 
line 312 (see Fig. 3b), which may be defined as 

30 cr = (27cB C RMsd)" 1 , (8) 



5 



10 
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in which "B C rms" is the centralized root mean square 
bandwidth of the transmitted pulse. The RMSE values 
continue to increase at a faster rate within region III (see 
Fig. 3b) . It is noted that B C rms may be expressed as 

5 # 

Bcrms = (J(f-fc) 2 P SD (f)df) 1/2 / (9) 

in which w fcf is the center frequency of the pulse, and the 
integration operation is performed from 0 to +oo. It is 
10 further noted that fc may be expressed as 

f c - JfP 8D (f)df, (10) 

in which the integration operation is performed from 0 to 
15 +oo. Moreover, the root mean square bandwidth may be 
expressed as 

Brms 2 = Bcrms 2 + f c 2 - (11) 

20 Accordingly , in the event the center frequency fc is much 
larger than the centralized root mean square bandwidth B C rms, 

Brms « fc. (12) 

25 The behavior of the simulation curves 301-307 within 

region I (see Figs. 3a-3b) is characteristic of . the 
performance of a "coherent" receiver, which estimates the 
range of a target relative to a peak of the ambiguity 
function within the width of the function's main lobe. The 

30 behavior of the curves 301-307 after their sharp transitions 
from lower RMSE levels to higher RMSE levels within region 
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II (see Figs. 3a-3b) is characteristic of the performance of 
a "semi-coherent" receiver, which estimates target range 
relative to the peak of the . envelope of the ambiguity 
function. .As illustrated in Fig. 3a, range estimates 
5 provided by the semi-coherent receiver have associated 
errors (RMSE) thatarfe significantly higher. than the errors 
associated with the range estimates of the coherent 
receiver. 

The embodiments disclosed herein will be better 
10 understood with reference to the following illustrative 
examples and Fig. 1. In a first example, the sonar 
transmitter 103 transmits a single sonar pulse ("ping" ) 
having a frequency range of about 87 kHz through a 
transmission medium such as water. The transmitted ping 
15 travels through the water until it strikes the selected 
target 112, which returns a reflected sonar pulse ("echo") 
to the sonar sensor 102. It is understood that a typical 
target may return multiple echoes. The target 112 in this 
example is described as returning a single echo for clarity 
20 of illustration. The sensor 102 provides a signal 

representative of the echo to the sonar receiver 104, which 
subsequently provides an indication of the echo to the 
plurality of band-pass filters 106. 

In this illustrative example, the plurality of band- 
25 pass filters 106 passes respective frequency sub-bands 
ranging from 6-18 kHz, 18-31 kHz, 31-43 kHz, 43-56 kHz, 56- 
68 kHz, 68-81 kHz, and 81-93 kHz, thereby spanning the ping 
frequency range of 87 kHz. The center frequencies fc of the 
frequency sub-bands are therefore approximately equal to 12 
30 kHz, 25 kHz, 37 kHz, 50 kHz, 62 kHz, 75 kHz, and 87 kHz, 
respectively. Further, the centralized root mean square 
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bandwidth B C rms of the transmitted pulse is fixed at about 
2.1 kHz. 

Next, the echo and a representation of the transmitted 
ping are filtered by each of the band-pass filters 106.1- 
5 106. n, which provide filtered versions of the echo, and ping 
* to the "cross correlator 108. The cross correlator 108 then 
performs multiple cross correlation operations on the 
. filtered -versions of the echo and ping, and provides 
multiple sets of cross correlation output data corresponding 
10 to the respective frequency sub-bands, to the data analyzer 

110., • 

As described above", the data analyzer 110 of the 
presently disclosed embodiment is operative to analyze the 
cross correlation output data to determine the variability 

15 of cross correlation peaks within each frequency sub-band, 
and to identify the lowest frequency sub-band having a 
corresponding low peak ambiguity. To this end, the data 
analyzer 110 effectively constructs an ambiguity function 
like the ambiguity functions 201-203 (see Figs. 2a-2b) for 

20 each of the frequency sub-bands 6-18 kHz, 18-31 kHz, 31-43 
kHz, 43-56 kHz, 56-68 kHz, 68-81 kHz, and 81-93 kHz, and 
determines the peak variability for each frequency sub-band 
using the ambiguity functions. 

In this illustrative example, it is assumed that an 

25 analysis of the cross correlation output data by the data 
analyzer 110 determines that the peak variability 
corresponding to each of the frequency sub-bands 56-68 kHz, 
68-81 kHz, and 81-93 kHz is low (i.e., the peak variability 
is within the width of the main lobe in each of the 

30 corresponding ambiguity functions) , and the peak variability 
corresponding to each of the frequency sub-bands 6-18 kHz, 
18-31 kHz, 31-43 kHz, 43-56 kHz is high (i.e., the peak 
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variability varies between the side lobes in each of the 
associated ambiguity functions) . In this example, the 
operation of the sonar system 100 is therefore like that of 
a coherent receiver when transmitting pings having 
5 frequencies ranging from 56-93 kHz, and a semi-coherent 
receiver when transmitting pings having frequencies ranging 
from 6-56* kHz. 

Accordingly, the data analyzer 110 identifies the 
frequency sub-band 56-68 kHz. as the lowest frequency sub- 

10 band having a corresponding low peak ambiguity. Next, the 
data analyzer "110 analyzes the peak variability data versus 
SNR for the identified sub-band 56-68 kHz, as depicted in 
Fig. 3a by the curve 305, and determines the break point on 
the curve where the sharp transition from low RMSE levels to 

15 high RMSE levels begins, i.e., break point 7. It is noted 
that break point 7 on the peak variability versus SNR curve 
305 corresponds to a RMSE of about 10~ 7 and an SNR of about 
30 dB. The data analyzer 110 therefore estimates the SNR 
within the sonar environment to be about 30 dB for echoes 

20 and pings with a center frequency of 62 kHz and a 
centralized root mean square bandwidth of 2.1 kHz. 

As described above, the data analyzer 110 provides 
information to the controller 111, which uses this 
information to cause the sonar transmitter 103 to modify the 

25- characteristics of the transmitted pings to adapt to the 
estimated SNR. To this end, the data analyzer 110 
determines that the peak variability versus SNR curve 305 
corresponds to the center frequency 62 kHz for the frequency 
sub-band 56-68 kHz. In this illustrative example, the 

30 center frequency 62 kHz is the highest center frequency that 
has not yet been affected by ambiguity. The data analyzer 
110 then identifies this frequency to the controller 111, 
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and the controller 111 causes the sonar transmitter 103 to 
emit pulses toward the target with a center frequency of 62 
kHz and above, up to the maximum frequency available to the 
system hardware. In this way, the sonar system 100 
5 effectively removes the frequency bands that would create 
large ambiguities. The - sonar system 100 then uses these 
pulses to obtain more accurate sonar range estimations 
within the environment of interest. For example, the system 
100 may determine, with low ambiguity, the location of the 
10 cross correlation peak' corresponding to the range of the 
selected target 112, and may then determine the range of the 
target 112 based on the peak location. 

It should be understock that instead of transmitting a 
single ping having a frequency range of about 87 kHz and 
15 band-pass filtering the resulting echo, as in the first 
example above, a single ping corresponding to each of the 
frequency sub-bands 6-18 kHz, 18-31 kHz, 31-43 kHz, 43-56 
kHz, 56-68 kHz, 68-81 kHz, and 81-93 kHz may alternatively 
be transmitted. Next, multiple cross correlation operations 
20 may be performed on the echoes and pings corresponding to 
each frequency sub-band, and multiple sets of cross 
correlation output data corresponding to the respective 
frequency sub-bands may be provided to the data analyzer, 
which generates results essentially the same as those 
25 generated in the first example above, i.e., . a single 
distribution is formed by one estimate per frequency band. 
In both cases, the center frequency fc may be deteinuined 
such that the overall distribution of cross correlation 
output data is the combination of two distributions having 
30 different standard deviations, i.e., a first distribution 
corresponding to estimates above the center frequency and a 
second distribution corresponding to estimates below the 
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center frequency, in which the second distribution has a 
standard deviation smaller than the first distribution. The 
determined center frequency fc corresponds to the breakpoint 
of interest, e.g., break point 7 in the first example above 
5 (see Fig. 3a) . It is also understood that the data analysis 
results. obtained from transmitting the single broadband ping 
and those obtained from transmitting the multiple narrow- 
band pings may be combined to determine the center frequency 
f c corresponding to the breakpoint of interest . 

10 . In a second example/ instead of transmitting a single 

broadband ping having a frequency range of about 87 kHz, the 
sonar transmitter 103 transmits multiple broadband pings 
through the water. It is understood that the centralized 
root mean square bandwidth B C rms of the multiple pings is 

15 fixed at about 2.1 kHz. As in the first example described 
above, the multiple pings travel through the water until 
they strike one or more selected targets 112, which return 
reflected sonar pulses (echoes) to the sonar sensor 102. 
The sensor 102 in turn provides signals representative of 

20 the echoes to the sonar receiver 104, which subsequently 
provides indications of the echoes to the plurality of band- 
pass filters 106. 

As in the first example above, the plurality of band- 
pass filters 106. 1-106. n passes respective frequency sub- 

25 bands ranging from 6-18 kHz, 18-31 kHz, 31-43 kHz, 43-56 
kHz, 56-68 kHz, 68-81 kHz, and 81-93 kHz. The center 
frequencies of the respective frequency sub-bands are 
therefore approximately equal to 12 kHz, 25 kHz, 37 kHz, 50 
kHz, 62 kHz, 75 kHz, and 87 kHz. Accordingly, the sonar 

30 receiver 104 provides the pings having center frequencies 12 
kHz, 25 kHz, 37 kHz, 50 kHz, 62 kHz, 75 kHz, and 87 kHz to 
the band-pass filters 106. 1-106. n passing sub-bands 6-18 
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kHz, 18-31 kHz, 31-43 kHz, 43-56 kHz, 56-68 kHz, 68-81 kHz, 
and 81-93 kHz, respectively. 

The echoes and representations of the corresponding 
transmitted pings are then filtered by the band-pass filters 
5 106, which provide filtered versions of the echoes and pings 
to the cross correlator 108. The cross correlator 108 
performs ^multiple cross correlation operations on the 
filtered versions of the echoes and' pings, and provides 
multiple se.ts of * cross correlation output data to the data 
10 analyzer 110. 

Next, the data analyzer 110 effectively constructs an 
ambiguity function for the~ echoes and corresponding 
transmitted pings. For example, Fig. 4a depicts 

representative ambiguity functions for multiple broadband 
15 pings corresponding to frequency sub-band 31-43 kHz, Fig. 4b 
depicts representative ambiguity functions for multiple 
broadband pings corresponding to frequency sub-band 56-68 
kHz, and Fig. 4c depicts representative ambiguity functions 
for multiple broadband pings corresponding to frequency sub- 
20 band 81-93 kHz. In each of the ambiguity functions depicted 
in Figs. 4a-4c, a black dot is employed to indicate an 
estimated location of a cross correlation peak, as computed 
by the data analyzer 110. Accordingly, the data analyzer 
110 may use the estimated peak locations (as indicated by 
25 the black dots) to perform a statistical analysis for 
determining the level of peak ambiguity for each ambiguity 
function. 

As shown in Figs. 4a-4c, the peak variability for each 
one of the pings having a center frequency equal to 37 kHz, 
30 62 kHz, and 87 kHz, respectively, is indicated by the 
estimated locations of the corresponding cross correlation 
peaks. For example, for the pings corresponding to 
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frequency sub-band 81-93 kHz (fc = 87 kHz), the majority of 
the computed peaks are located within the width of the main 
lobe of the ambiguity function when SNR = 23 dB, thereby 
indicating low peak ambiguity. However, when the SNR is 
5 successively reduced from 23 dB to 18 dB, 13 dB, and 8 dB, 
the peak ambiguity gradually increases, as indicated by the 
increasing number of computed peaks located outside the main 
lobe (see Fig. 4c) . Similarly, for the pings corresponding 
to frequency sub-band 56-68 kHz (fc « 62 kHz), there is low 
10 peak ambiguity when SNR « 23 dB, and increasing peak** 
ambiguity as the SNR is successively reduced from 23 dB to 
18 dB, 13 dB, and 8 dB (see Fig. 4b); and, for the pings 
corresponding to frequency sub-band 31-43 kHz (fc = 37 kHz), 
there is low peak ambiguity when SNR = 23 dB and 18 dB, and 
15 increasing peak ambiguity when the SNR is successively 
reduced from 18 dB to 13 dB and 8 dB (see Fig. 4c) . 

Accordingly, the technique of this second example may 
be employed to obtain an empirical estimation of the SNR in 
the environment of interest. For example, as shown in Figs. 
20 4a-4c, the lowest possible SNR is empirically estimated to 
be about 18 dB. It is understood that varying levels of 
peak ambiguity and SNR may be obtained depending on the 
center frequency and frequency range of the pings, and the 
specific characteristics of the sonar environment. The data 
25 analyzer 110 therefore estimates the SNR within the sonar 
environment to be about 18 dB for echoes and pings with a 
center frequency of 37 kHz and a centralized root mean 
square bandwidth of 2.1 kHz. 

Next, the data analyzer 110 provides information to the 
30 controller 111, which uses this information to cause the 
sonar transmitter 103 to modify the characteristics of the 
transmitted pings to adapt to the estimated SNR. In this 
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second example , the center frequency 37 kHz is the highest 
center frequency that has not yet been affected by 
ambiguity- The data analyzer 110 then identifies this 
frequency to the controller 111, and the controller 111 
5 causes the sonar transmitter 103 to emit pings toward the 
targets 112 with a center frequency equal to 37 kHz and 
above , up to the maximum frequency available to the system 
hardware. In this way, the sonar System 100 effectively 
removes the frequency bands that would create ^ large 
10 ambiguities. The sonar system 100 then uses these pings to 
obtain more accurate sonar range estimations within the 
environment of interest. ' For example, the system 100 may 
determine, with low ambiguity, the locations of the . cross 
correlation peaks corresponding to the range of the selected 
15 targets 112, and may then determine the respective target 
ranges based on the peak locations. 

It should be understood that instead of transmitting 
multiple broadband pings and band-pass filtering the 
resulting echoes, as in the second example above, multiple 
20 narrow-band pings corresponding to each of the frequency 
sub-bands 6-18 kHz, 18-31 kHz, 31-43 kHz, 43-56 kHz, 56-68 
kHz, 68-81 kHz, and 81-93 kHz may alternatively be 
transmitted, thereby obviating the need for band-pass 
filtering. In both cases, the results generated by the data 
25 analyzer are essentially the same, i.e., many estimates for 
each frequency sub-band, and therefore many distributions 
for the peak estimates. It is understood that the data 
analysis results obtained from transmitting the multiple 
broadband pings and the multiple narrow-band pings may be 
30 combined for empirically estimating the lowest possible SNR. 

A first method of operating the presently disclosed 
system for estimating the SNR in a noisy sonar environment 
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is illustrated by reference to Fig. 5a. As depicted in step 
502 , a sonar transmitter transmits a single sonar pulse 
("ping") having a predetermined frequency range through a 
transmission medium. Next, a sonar sensor/receiver 

5 . receives, as. depicted in step 504, a reflected sonar signal 
(" echo" ) reflected from a selected target. The receiver 
then provides, as depicted in step 506, an indication of the 
echo to a plurality of band-pass filters, which pass 
respective frequency sub-bands substantially spanning the 

10 predetermined frequency range. Next, the band-pass filters 
filter, as ■ depicted in step 508, the echo and a 
representation of the transmitted ping. A cross correlator 
then performs multiple cross correlation operations on the 
filtered versions of the echo and pulse, as depicted in step 

15 510, and provides multiple sets of cross correlation output 
data corresponding to the frequency sub-bands to a data 
analyzer. Next, the data analyzer determines, as depicted 
in step 512, the variability of cross correlation peaks 
within each one of the frequency sub-bands, and identifies, 

20 as depicted in step 514, the lowest frequency sub-band 
having a corresponding low peak ambiguity. The data 
analyzer then determines, as depicted in step 516, the break 
point on the peak variability versus SNR curve corresponding 
to the identified frequency sub-band in this sonar 

25 environment where the sharp transition from low error levels 
to high error levels begins. Next, the data analyzer 
estimates, as depicted in step 518, the SNR in the sonar 
environment to be the SNR corresponding to the break point 
determined in step 516. The data analyzer then determines, 

30 as depicted in step 520, the center frequency corresponding 
to the identified frequency sub-band. This is the highest 
center frequency that has not yet been affected by 
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ambiguity. Next, the data analyzer identifies the center 
frequency to a controller, as depicted in step 522, and the 
controller controls the sonar transmitter to emit pulses 
toward the target at frequencies greater than or equal to 
5 that center frequency, up. to the maximum frequency available 
to the system hardware. - The sonar system then uses these • 
pulses, as depicted in step 523, to obtain more accurate 
sonar range estimations within the environment of interest. 

A second method of operating the presently disclosed 
10 system for estimating the SNR in a noisy sonar environment 
is illustrated by reference to Fig. 5b. As depicted in step 
524, a sonar transmitter transmits multiple broadband pings 
having a predetermined frequency range and predetermined 
center frequencies fc through a transmission medium. Next, 
15 a sonar sensor/receiver receives, as depicted in step 525, 
reflected sonar signals (echoes) reflecting from one or more 
selected targets. The receiver then provides, as depicted 
in step 526, indications of the echoes to a plurality of 
band-pass filters, which pass respective frequency sub-bands 
20 having respective center frequencies. Next, the band-pass 
filters filter, as depicted in step 528, the echoes and 
representations of the transmitted pings. A cross 

correlator then performs multiple cross correlation 
operations on the filtered versions of the echoes and pings, 
25 as depicted in step 530, and provides multiple sets of cross 
correlation output data corresponding to the respective 
center frequencies to a data analyzer. Next, the data 
analyzer computes, as depicted in step 532, a cross 
correlation peak for each ping, and determines, as depicted 
30 in step 534, the variability of cross correlation peaks 
within each sub-band by a statistical analysis of the peak 
locations. The data analyzer then determines, as depicted 
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in step 536, the center frequency that allows the lowest SNR 
while maintaining a low peak ambiguity, using the results of 
the statistical analysis of step 534. This is the highest 
center frequency that has not yet been affected by 
5 ambiguity. Next, the data analyzer identifies the center 
frequency to a controller, as depicted in step 538, and the 
controller controls the sonar . transmitter to emit pulses 
toward the target at frequencies greater than or equal to 
that . center frequency, up to the maximum frequency available 
10 to the system, hardware. The sonar system then uses these 
pulses, as depicted in step 540, to obtain more accurate 
sonar range estimations within the environment of interest. 

Having described the above illustrative embodiments, 
other alternative embodiments or variations may ^ be made.. 
15 For example, although the controller 111 is depicted in Fig. 
1 as a separate functional block, it should be recognized 
that the functionality of the controller 111 may 
alternatively be integrated either within the data analyzer 
110 or within the transmitter 103. 
20 It should also be understood that when the 

characteristics of the transmitted pulse are unknown, e.g., 
in the case of a passive sonar system, the transmitted pulse 
may be estimated by averaging using multiple echoes or 
returns. Such averaging may take into account some 
25 distortions (known or estimated) to the pulse. The 
estimated signal may then be applied to the cross 
correlator, as described above. 

In addition, it was described in the second example 
above that multiple echoes and pings are provided to a 
30 plurality of band-pass filters for filtering before 
performing cross correlation operations on the echoes and 
pings. However, in an alternative embodiment, the plurality 
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of filters may be omitted, and unfiltered representations of 
the echoes and pings may be provided directly to the cross 
correlator. 

In addition, it was described that for .pulses having a . 
5 fixed centralized root mean square bandwidth B C rms and 
different* center frequencies fc, the behavior, of the 
simulation curves 301-307 witfein region I (see Figs. 3a-3b) 
is characteristic of the performance of a coherent receiver, 
and the behavior of the curves 301-307 following their sharp 

10 transitions from low error levels to high error levels 
within region II (see Figs. 3a-3b) is characteristic of the 
performance of a semi-coherent receiver. However, it should 
be understood that simulation curves similar to the curves 
301-307 may be generated to determine the root mean square 

15 error (RMSE, seconds) as a function of SNR (dB) and 
centralized root mean square bandwidth B C rms for a fixed 
center frequency fc. 

In this case, the behavior of the simulation curves 
would be characteristic of the performance of a coherent 

20 receiver for high SNR levels, and as the SNR level 
decreases, the pulses with lower B C bms values would be 
affected by peak ambiguity first, thereby causing their 
performance to degrade to that of a semi-coherent receiver. 
In general, the pulses with the larger B C rms values are more 

25 resilient to peak ambiguity. As shown in Figs. 3a-3b, as 
SNR levels decrease, the pulses with higher center 
frequencies fc are affected by peak ambiguity first, and the 
pulses with lower center frequencies fc are more resilient 
to peak ambiguity. In a practical sonar system, however, 

30 only a limited range of frequencies may be employed for the 
sonar pulses and the sonar receiver. This limitation places 
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bounds on the pulses' center frequency fc and centralized 
root mean square bandwidth B C rms* 

The effect of this limited frequency range within the 
practical sonar system is illustrated by the performance 
5 curve 308 (see Fig. 3b) . For example, the maximum allowable 
Bcbms (e.g., 2.1 kHz) is the same for all of the pulses used 
to generate the curve 308. Further, for high levels of SNR 
(region I; see Fig. 3b), the behavior of the performance 
curve 308 is characteristic of a coherent receiver up to 
10 break point 9, which corresponds to the maximum allowable "fc* 
(e.g., 87 kHz). These maximum allowable B C rms and fc values 
are a direct consequence of the limited frequency range of 
the practical sonar system. 

In addition, the sonar system 100 (see Fig. 1) may 
15 further include a user input device such as a keyboard or 
control panel to allow the user to configure the system. 
Further, the sonar sensor 102 may include at least one 
suitable sonar transducer (e.g., a hydrophone sensor) 
operative to detect an echo signal, and the sonar receiver 
20 104 may include conventional filters and amplifiers for 
enhancing the echo signal before providing the signal to the 
plurality of band-pass filters 106. The sonar system 100 
may also include one or more computers operative to store 
and process the sonar signal data and the cross correlation 
25 output data. For example, the data analyzer 110 and the 
controller 111 may include one or more microprocessors, 
application specific integrated circuits (ASICs) , and/or 
microcomputers operative to analyze the signal and output 
. data and to control the sonar transmitter 103, in accordance 
30 with the methods disclosed herein. In addition, it is 
appreciated that the presently disclosed system for 
estimating SNR in noisy environments may be implemented for 



WO 2005/006009 PCT/US2004/0 18059 

simulation and/or operation using one or more programmed 
general purpose computers and/or special purpose hardware. 

It should be understood that the presently disclosed 
systems and methods are suitable for measuring delays 
5 between a set of consecutive signals , such as echoes or 
returns from a complex target or multiple targets. These 
signals may be pulses provided by an active transmission 
source such ' as an active sonar system/ or may comprise a 
transmitted pulse and an echo returning from the selected 
10 target. 

It will also be appreciated by those of ordinary . skill 
in the art that modifications to ' and variations of the 
above-described noise adaptive sonar signal processor may be 
made without departing from the inventive concepts disclosed 
15 herein. Accordingly, the invention should not be viewed as 
limited except as by the scope and spirit of the appended 
claims . 
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